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ABSTRACT 

We study the general evolution of spherical over-densities for thawing class of dark 
energy models. We model dark energy with scalar fields having canonical as well as 
non-canonical kinetic energy. For non-canonical case, we consider models where the 
kinetic energy is of the Born-Infeld Form. We study various potentials like linear, 
inverse-square, exponential as well as PNGB-type. We also consider the case when 
dark energy is homogeneous as well as the case when it is inhomogeneous and virializes 
together with matter. Our study shows that models with linear potential in particular 
with Born-Infeld type kinetic term can have significant deviation from the ACDM 
model in terms of density contrast at the time of virialization. Although our approach 
is a simplified one to study the nonlinear evolution of matter overdensities inside the 
cluster and is not applicable to actual physical situation, it gives some interesting 
insights into the nonlinear clustering of matter in the presence of thawing class of 
dark energy models. 
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1 INTRODUCTION 

One of the most significant discoveries in cosmology in re- 
cent years is the fact that our universe is currently going 
throu gh an accelerated expansion phase (|Riessll2004l : iKnop I 
I2003D . This can have far reaching implications for funda- 
mental theories of physics. This late time acceleration of 
the universe can be due to the presence of an exotic fluid 
with large negative pressure known as dark energy or due 
to the modification of gravity itself. The simplest candi- 
date of dark energy is provided by cosmological constant 
A with equation of state parameter w = —1. However, 
the ACDM model is plagu ed with fi n e tun in g and cosmi c 
coincidence problems (See ICopelandl |2006l); ISamJ ||2009|): 
ISahni fc Starobinskv I feoOOl: Padmanabhan ( 2003) : iLindej 
l|2008l ): lFriemanl (|2008l ): iCaldwelli ()200a ): iSilvestrii (\200^ for 



a nice review j. 

Scalar field models mimicking a variable A can 
alleviate the fine tuning and coincidence problems 
and provide an interesting al terna tiv e to co s molog- 
ical constant (|Ratra fc Peebl3 Il988l : ICaldwelll 1 19981 : 
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iLiddle fc Scherrej 1 19991 : ISteinhardd Il999l l. These scalar 
field models can be broadly classified into two categories 
depending upon the form of their potentials: fast roll and 
slow roll mod els termed as freezi ng and thawing models in 
the literature ()Caldwell fc Lindeil [2005). In case of fast roll 
models, the potential is steep resulting the scalar field to 
track the background fiuid and is sub dominant for most of 
the evolution history. Only at late times, the field becomes 
dominant and drives the acceleration of the universe. Such 
solutions are also known as trackers. 

Then there are the slow-roll models for which the field 
kinetic energy is much smaller than its potential energy. 
Usually it has sufficiently flat potential similar to an inflaton. 
At early times, the field is nearly frozen at w = — 1 due to the 
large Hubble friction. Its energy density is nearly constant 
with a negligible contribution to the total energy density of 
the universe. But as radiation/matter rapidly dilutes due to 
the expansion of the universe and the background energy 
density becomes comparable to scalar field energy density, 
the field breaks away from its frozen state and evolves slowly 
to the region with w > —1. However, in this case, the model 
needs some degree of fine tuning of the initial conditions in 
order to achieve a viable late time evolution. 

Recent observations suggest that the equation of state 
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parameter for dark energy does not signific antly deviate 
from w = — 1 a round the present epoch (|Wood-Vasevl 
I2OO7I : lDavisll2007l ). This type of equation of state can be 
easily obtained in dynamical models represented by thaw- 
ins scalar fields. Motiv ated by this fact, Scherrer and Sen 
(|Scherrer k, Senll2008al ) examined quintessence models with 
nearly flat potentials satisfying the slow-roll conditions. It 
was shown that under the slow-roll conditions, a scalar field 
with a variety of potentials V{(j)) evolves in a similar fash- 
ion and one can derive a generic expression for the equa- 
tion of state for all such scalar fields. Simi lar results were 
later established for the case of phantom JScherrer fc SenI 
l2008bl ') and tachyon scalar fields (|Aiill2009l : [Senll2010l '). It was 
demonstrated that under slow-roll conditions, all of them 
have identical equation of state and hence can not be distin- 
guished, atleast, at the level of background cosmology. The 
crucial assumption for arriving at this important conclusion 
was the fulfillment of the slow-roll conditions for the field 
potentials. 

In a recent paper by Sen etal. l|Senl2010l '). assumption of 
slow roll has been relaxed, but it was assumed that the scalar 
field is of thawing type i.e it is initially frozen at ui = — 1 
due to large Hubble damping. With these choices the evolu- 
tion of a variety of scalar field models having both canonical 
and non-canonical kinetic terms have been studied. Obser- 
vational quantities like Hubble parameter, luminosity dis- 
tance and quantities related to the Baryon Acoustic Oscil- 
lation (BAO) measurement have been studied for variety of 
potentials. It was shown that given the current error bars 
for different observational data related with background cos- 
mology, it is difficult to distinguish different thawing models 
from ACDM. It was also shown that lower the value of Q.mO 
parameter, higher the chance to distinguish and there is a 
redshift range between z ~ 0.4— 1 where the deviation from 
ACDM is maximum irrespective of the model considered. 

Dark energy not only affects background expansion 
rate and thereby modifying the observables like distance- 
redshift relation, it also affects the growth of struc- 
ture. On large scales, dark energy also clusters which in 
turn affects the matter clustering. But at small scales, 
although dark energy is more or less smoothly dis- 
tributed, it provides an extra Hubble drag thereby slow- 
ing down the growth rate of matter perturbations. The 
behavior of linear perturbation in a scalar field and 
its effect on large scale structure of th e univer s e hav e 
been studied by a number of author s llBartold ( 2004 ): 
Hul (l2005h; iGordon fc Wandj (120051} : lAvelind l2008li: 



Unnikrishnaml (l2008a^ : 



■ ■ .Gordon fc Hul J2OO4): IZimdahll 

( 20051 ): lUnnikrishnanI ((2008bl ): |Jassa'iT(|2009l . I2OI0I )). The 
impact has been studied for cosmic microwave back- 



ground (M acor r;^ 20031: Baccigalupi fc Acguavival I2OO6I : 
lAcauaviTO fc Baccigalu pi 2006: iGiovil I200"), gala x y red - 
shift surveys jEisenstein (. loos') Hseo fc Eisensteiiil l2003h: 
iHaiman I (|2000| ): IJairi fc Taylor I (|2003l '): iBernstein fc JainI 
(|2004l ): lAllenI (|2004l ')). cross-correlatio n of the inte- 
grated Sachs- Wolf effect (|Pogosian 2005: ICorasanitil I2OO5I : 
iGiannant onio 2008has well as in neutral hydrogen surveys 
(|Bharadwai..2009, '). 

In the non-linear regime one has to perform N-Body 



simulations to accurately measure the effect of dark energy 
on the clustering properties of ma tter and this has been in- 
vestigated in a number of papers (|Maccidl200j : lBaldill2008l : 
Courtin 20IQ). The simplest analytical approach to study 
the nonlinear clustering of matter is the sph erical collapse 
form alism first developed by Gunn and Gott (|Gunn fc Gotd 
Il972h . It describes how a small spherical over density 
decouples from the background evolution, slows down, then 
eventually turns around and collapses. It is generally as- 
sumed that the collapse is not complete so the system does 
not reach the singularity, instead it eventually virializes 
and stabilizes to a finite s ize. Combining this wit h the 
Press- Schechter formalism (jPress fc Schechteil 1 19741 '). one 
can have a model for formation of structures which predicts 
the abundances of virialized objects as a function of mass. 
This model can give insight into the physics of structure for- 
mation which can be ultimately used for a detailed N-Body 
simulation. There are numerous studies of the spherical 
collapse formalism to inclu de both homogeneous and 
inhomogeneous dark en e rgy ( Lahavl I1991I: IShaw fc Motal 
I2OO8I; I Manera & Mot4 
Basilakos fc Voglis 2 0071: 
B artelmannI bood: 



200e : iLokas & Hoffman 



2001; 
2004- 



- - - , . Schacfo r fc Koyama _ . 

Basilakoj I2OO3I: iHorellou fc Berge I 



20051: iMaor fc Lahavl l2005l: iMota fc Van de Bruckl | 20ol : 



Nunes fc Motall200d: IWang fc Steinhardtlll998l : IWandboOQ : 
Weinberg fc Kamionkows ki"2003':'Zeng| |2005l ') 



In this work, we systematically study the spherical col- 
lapse model for thawing type scalar field dark energy models 
for various potentials. We consider the canonical as well as 
noncanonical kinetic energy for the scalar fields. We assume 
the dark energy to be homogeneous as well as inhomoge- 
neous when it also virializes together with matter inside the 
spherical over density. We calculate the matter density con- 
trast at turnaround as well as at virialization and compare 
the result with the corresponding ACDM model. 

The plan of the paper is as follows: in section 2 we 
describe the background evolution of thawing dark energy 
models for ordinary scalar field as well as for tachyons. In 
section 3 we study the evolution of linear and nonlinear mat- 
ter density perturbation assuming the dark energy to be ho- 
mogeneous. In section 4 we describe the spherical collapse 
formalism for a general scalar field dark energy model and 
study the virialization process in this scenario and discuss 
the results obtained. Finally, we draw our conclusions in 
section 5. 



2 BACKGROUND EVOLUTION FOR 
THAWING DARK ENERGY MODELS 

As discussed in the introduction, thawing dark energy mod- 
els are characterizedin such a way that in the early universe 
the scalar field is frozen by very large Hubble damping due 
to the expansion of the universe. As the universe expands the 
Hubble parameter decreases so the Hubble damping and the 
scalar field starts evolving slowly down its potential. There- 
fore, the equation of state initially starts with w = —1 and 
slowly departs from this value in the later time. In this paper 
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we consider both ordinary scalar field with canonical kinetic 
term as well as tachyon type scalar field having Born-Infeld 
type kinetic term. 



2.1 Thawing Scalar field 

In what follows, we shall assume that the dark energy is de- 
scribed by a minimally-coupled scalar field, 0, with equation 
of motion 



if) + 3H(j} + dV/d(j) = 

where the Hubble parameter H is given by 



H = 



(1) 



(2) 



Here p is the total energy density in the universe. We model 
a flat universe containing only matter and a scalar fleld, so 
that il^ + Q,M = 1. 

Equation U]) indicates that the fleld rolls downhill in 
the potential V((j}) but its motion is damped by a term pro- 
portional to H. The equation of state parameter w is given 
hy w — / where the pressure and density of the scalar 
field have the forms 



P0 = y-^^W- 



P4> 



(3) 
(4) 



Defining A 



\ 2 



one can 



form an autonomous system of equations involving the two 
observables fl^ and 7 = (l + w) to gether with the parameter 
A. This is derived by jScherrer fc Sen,2008a )as 



7 



: -37(2-7) + A(2-7)y3^, 
: 3(l--7)f7^(l-n^), 

A' = -^A^(r-i)y^. 



(5) 
(6) 
(7) 



Given the initial conditions for 7, fi^ and A, one can 
solve this system of equations numerically for different po- 
tentials. We are interested in thawing models i.e models for 
which the equation of state is initially frozen at w = —1. 
Hence, we take 7 = initially for our purpose. Also one has 
to choose the initial value Ai. It was earlier shown that if Xi 
is small (which is equivalent to assuming the slow-roll condi- 
tions similar to inflation) the scalar fleld evolution is simila r 
to that of cosmological constant A (jScherrer fc Senll2008al ). 
Here we assume that the slow-roll conditions are strongly 
broken i.e Ai = 1 in order to have a behavior which is in 
principle different from A. In general the contribution of the 
scalar field to the total energy density of the universe is in- 
significant at early times, nevertheless one has to fine-tune 
the initial value of £7^ in order to have its correct contribu- 
tion at present. This is the fine tuning one needs to have 
in a thawing model. With these initial conditions we evolve 
the above system of equations from redshift 2 = 1000 (or 
a — 10""^) till the present day z = (a = 1). We consider 
various types of potentials e.g V = 4>, V = 4>^ , V = e"^ and 
V = 4>~^ , characterized by F = 0, |, 1 and | respectively. 



We assume fJ^o ~ 0.75 at the present epoch for all chosen 
values of V. 

We also consider the Pscudo-Nambu Goldstone Boson 
(PN GB) model (jFrieman 1995). (For a recent discussion, see 
Ref. llAbrahamsdbOOTh and references therein). This model 
is characterized by the potential 

l/(0) = m^[cos(<^//) + l]. (8) 

Alam et al., l|Alamll2003l ) have previously considered such 
type of potential to see whether dark energy is decaying or 
not. 

Without any lose of generality, we choose f = 1. Here, 
the constant m is related to the mass of the fleld. 



2.2 Thawing Tachyon Model 

In recent years, tachyon-like scalar field having ki- 
netic energy of Born-Infeld form has generated lot 
of interests in cosmology. There have been several 
investigations usi ng t hi s type of field as dark en- 



ergy (ICopelandl l2005l: pTsuiikawa fc Samil |2004|: ISam 



200 i iFeldeij |2002(: lAbramo fc Finellil l2003l : iBaela I 12003 
lAguirreeabiria fc Lazkozlliooi ). In what follows, we form 
the similar set of equations for tachyon field as described 
above for ordinary scalar field. 

The tachyon fiel d is specifi ed by the Dirac-Born-Infeld 
(DBI) type of actio nllSenll2002allbl: [Garousill2000l : iBergshoefj 
I2OOOI : iKluson |[2OO0l : iKutasov fc Niarchos 11200^ ): 



(9) 



In FRW background, the pressure and energy density of the 
tachyon field cf> are given by 



P<t> = 



-v(<^)yr 



The equation of motion which follows from ((9]) is 

') = 



4> + iH^{l~^') + ^{l 



(10) 
(11) 

(12) 



where H is the Hubble parameter and here prime denotes 
the derivative with respect to (j) ■ Now assuming A 



dV 



and F = V 



d4:''^ 
ldV_\2 

^ d(b ' 



we can also form an autonomous system 



of equatio ns involving O,^, 7 = (l + w) and A for tachyon 
(|Alill2009l '): 



7' = -67(1 - 7) + 2v/3^A(l - 7)3 
J7; = 3f7^(l-7)(l-n^) 



A' = -y3^A^(l-7)4(F--) 



(13) 
(14) 

(15) 



We choose the similar set of initial conditions to solve 
this system of equations for tachyon as we describe earlier 
for ordinary scalar field. We also assume the same set of 
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0.0 0.5 1.0 1.5 2.0 

Z 

Figure 1. Plot of equation of state w vs. redshift for different 
scalar field and tachyon models. Solid curves represent different 
Tachyon models with V{(j)) = (j),<t>^, e'^, (j>~^ respectively from top 
to bottom, Dashed curves from top to bottom represent different 
scalar field models with same potentials as in tachyon. Dotted 
curve represents PNGB model. HmO = 0.25. 

potentials V = (j), V = (f)^ , V = e'^ and V = (f>~^, charac- 
terized by r = 0, ^, 1 and | respectively and f^^o = 0.75. 



2.3 Background Result 

Let us now see the behavior of the different dark energy mod- 
els that we have considered above. In Figure 1, we plot the 
behaviors of the equation of state parameter w for different 
thawing models. It shows that the equation of state of differ- 
ent fields with different potentials behave differently as one 
approaches the present day although in the past their behav- 
iors are almost identical. This is not surprising as we have 
assumed the violation of slow-roll condition, i.e Xiniuai ~ 1- 
With slow-roll condition satisfied, i.e, Xi « 1, it was shown 
earlier that models with different potential s have the iden- 
tical w j a) both fo r scalar and tachyon fields (|Scherrer fc SenI 
l2008al lbl: lAiill2009l '). Although we show equations of state for 
some models which are more than w = —0.8 at present, 
these models are practically ruled out by current observa- 
tional data. 



2.0 F 




0.0 0.2 0.4 0.6 0.8 1.0 

a 



Figure 2. Variation of linear and nonlinear density contrast S 
as a function of scale factor for scalar field and tachyon models. 
Dashed line is for ACDM model, solid line is for scalar field with 
linear potential and dash-dotted line is tachyon field with linear 
potential. Dotted line is for scalar field with PNGB type potential. 
The upper set is the solution for full nonlinear equation given by 
eqn (16), where as the lower set is for linearized equation (17). 
C„0 = 0.25. 

One can solve these two equations numerically. To fix the ini- 
tial conditions, we assume that in early times,the universe is 
matter dominated with negligible dark energy contribution 
which is typical for thawing model. For the matter domi- 
nated regime, 5 ~ a and ^ ~ 1 fixes our initial conditions. 
In figure 2, we plot the evolution of both linear as well as 
nonlinear density contrast as a function of scale factor. We 
show the behaviors for scalar as well as tachyon field with 
linear and PNGB potential. To compare the results with 
ACDM, we also plot the corresponding behavior for ACDM 
case. As one can see from this plot, that although the lin- 
ear matter density contrasts for different thawing models do 
not deviate much from that of the ACDM model but the 
nonlinear density contrast has substantial deviation from 
ACDM model for both scalar as well as tachyon field. We 
should mention that although we show the results for linear 
and PNGB potentials, for other potentials, the behaviors 
are similar but with lesser deviations from ACDM model. 



3 EVOLUTION OF DENSITY 
PERTURBATIONS 

The matter density contrast is defined hy 5 — (pmc — 
Pm)/pm, where pmc is the perturbed matter energy density 
and pm is the background matter energy density. Assuming 
that the dark energy does not cluster, the evolution of the 
matter density contrast is governed by the equation 

5 + 2-5-47rGp™(l-h5)5- i— ^ =0. (16) 
a 3(1-1-0) 

In the linear regime, one can ignore higher order term in the 
above equation and approximate it as: 

5 + 2-5-4-kGp^5^0. (17) 
a 



4 THE SPHERICAL COLLAPSE MODEL 

The spherical collapse model, which has a long history in 
cosmology, is a simple and a fundamental tool for under- 
standing how a small spherical patch of over density forms 
a bou nd system via gravitational instability (|Gunn fc Gotj 
1 19721 ). As we consider flat, homogeneous and isotropic cos- 
mologies, driven by non relativistic matter and dark energy 
with the equation of state, — w{a)p^ with < 0, the 
equations that describe our background universe are given 
by 

SttG, 

[-) =— (Pm+P^) (18) 

and 
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Figure 4. Over density at turnaround vs. collapse redshift {zi-di) 
for different tachyon field models with the same potentials as in 
fig.3 for ClmO = 0.25. The lines correspond to the same as in fig.3. 



Figure 3. Over density at turnaround vs. collapse redshift (^coll) 
for different scalar field models with different potentials V{(t>) = 
4i,<j)'^,e'^,(j>~^ for Q^o = 0.25. In each figure, solid curve repre- 
sents inhomogeneous case, dashed curve represents homogeneous 
case and dotdashed represents ACDM model. Scalar field PNGB 
model is also compared with ACDM model in the bottom left 
figure. 



- = —[{3w{a) + l)p^ + pm] 



(19) 

where a{t) is the scale factor, pm = pmoaT^ is the back- 
ground matter density and p^ = P4,of{a) is the dark energy 
density, with 



/(o) = Exp[3 



1 + 



dw]. 



(20) 



The Priedmann equation can be written more simply as 
= {a/af = HQE(af and Ho is the Hubble constant 
with 

E{a) = [n^oo-^ + O^o/(a)] , (21) 

while f2mo = (87rGpmo)/3-ffo is the matter density param- 
eter and Sl^o = {8tvGp^o)/SHq is the corresponding dark 
energy parameter at the present epoch with Omo + fi(/>o = 1- 
The ^im{a) and f2,^(o) evolve with the scale factor as 

f2m(o) = and il^(a) = 



E^a) 



(22) 



Once we know the solution for f2^(a)and 7(a) by solving 
either (5)-(7) or (13)-(15), we can easily find the behavior 
of the normalized Hubble parameter in terms of fi,^ as 



E\a) 



HI 



(23) 



and the equation of state of dark energy w{a) = 7(a) — 1. 

The evolution of a spherical over dense patch of radius 
R(t) in the presence of dark energy is given by the Ray- 
chaudhuri equation: 



R 
R 



-AttG 



(w{R) 



1\ 1 



(24) 



where p^c and pmc are the dark energy density and the mat- 
ter density inside the spherical cluster respectively. w(R) 



represents the equation of state of the dark energy inside 
the cluster. With some initial over density, the spherical over 
dense region will expand until it reaches the maximum ra- 
dius (turn around point, R = 0) and then begins to collapse. 
After that the sphere virializes forming a bound system. 
With different initial density, it happens locally at different 
regions of the universe. Now, after performing the following 
transformations: 



a , R 
X = — and y = -=r 

at Rt 



(25) 



the equation of background evolution and that of the spher- 
ical perturbation become: 



X = Ht Q,m,t[^m{,x)x]~ 



and 

y = ■ 

where 

I{x,y) 
with 



^+i'yl{x,y) 



[l + 3wiR{y))]i^ 
[1 + 3w{x)] fix) 



(26) 



(27) 



Clustered DE 
Homogeneous DE 



(28) 



1-n, 



m,t 



Pm,t 



m,t 



(29) 



Subscript "t" denotes the turn around time. In order 
to solve the above set of equations, we have used the fact 
that the mass of the forming cluster is conserved: Pmc-R^ = 
Pmc,tRt and define the relation: 

i? \ CPm,t 



Pmc — Pmc,t 



(I)" 



(30) 



Here is the matter density contrast at turnaround which 

is defined as 

i^e^ = (^y\ (31) 



pm,t 



The function R{y) is given by R{y) = Rty = C ^^^o,ty and 
^m{x) by 



iiin{x) 



1 + vx'-''f{x) ' 



(32) 
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Figure 5. Plot of the non-linear density contrast at virialization 
as the function of collapse redshift (^con) fof different Scalar 
field models with V{4>) = 4>,<jP ,e'^ ,cf>~^ for = 0.25. In each 

figure, solid curve represents inhomogeneous case, dashed curve 
represents homogeneous case and dotdashed represents ACDM 
model. Scalar field PNGB model is also compared with ACDM 
model in the bottom last figure. 



Typically the scalar fields we consider here, have extremely 
small masses (of the order of present day Hubble scale in nat- 
ural units). This is necessary to have nearly flat potentials 
for the scalar fields at present day. Due to this, the scale of 
fluctuations for these scalar fields are extremely large, mak- 
ing it a smoothly distributed field within the horizon scale. 
So it is safe to assume the dark energy to be homogeneous. 
But it is still interesting to consider the case where the dark 
energy clusters along with the dark matter a nd avoid the en- 
ergy non- conservation problem examined in (|Maor fc Lahavl 
I2OO5I ). Hence in our subsequent calculations, we assume both 
the cases mentioned in equation (28). 

When dark energy is clustered together with the mat- 
ter, one can solve analytically the system of equations (26) 
and (27)as the function I(x,y) in eq.(28) depends only on 
y. So we easily calculate the density contrast C, at the turn 
around epoch by integrating the equations (26) and (27) 
and applying the boundary conditions {dy /dx)x^i ~ and 
y^^o ~ 0. The corresponding integral equation which gov- 
erns the behavior of ^, for a flat cosmological model is given 
below: 



1/2 



dy ■ 



(; + uyP{y) - {( + P{l)i^)y _ 

[x^lmix)]^''^ dx 



(33) 



In the case of homogeneous dark energy, one can not 
follow this analytical procedure and has to solve the system 
of equations (26) and (27) numerically with the boundary 
conditions {dy/dx)x=i ~ and yx=i ~ 1. 

In figure 3 and 4, we show the matter density con- 
trast at turnaround, ^ as a function of collapsed redshift 
for scalar and tachyon models with different potentials. To 
compare our result with ACDM model, we also show the cor- 
responding ^ for ACDM model. At higher Zcoii, overdensities 
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Figure 6. The plot of the non-linear density contrast at viri- 
alization as the function of collapse redshift(2:;,o;;) for different 
tachyon field models for QmO = 0.25 with the same potentials as 
in fig. 5. The lines correspond to the same as in fig. 5. 



at turnaround for different models tend toward the fiducial 
value of 5.6 for Einstein-de Sitter universe. For structures 
which are collapsing around present time i.e Zcoii ~ 0, the 
over densities at turnaround for different scalar and tachyon 
models are higher than the corresponding ACDM value. Also 
the values for the homogeneous case are slightly higher than 
the inhomogeneous case and the difference is more promi- 
nent for tachyon model with linear potential. For other cases 
the differences are extremely small. 

Just to mention, we show all the plots assuming flmo ~ 
0.25. For other values of Qmo, the behaviors are similar with 
lesser deviation from ACDM for higher fimo and vice versa. 

4.1 VIRIALIZATION IN THE SPHERICAL 
MODEL 

In general the spherical collapse formalism leads to a point 
singularity as the final state of the system. But physically 
the objects go through a virialization process and stabilize 
to a finite size. Such process of virialization is not built in 
the spherical coUpase formalism but we have to put it by 
hand in order to ensure virialization. 

There are number of things that enter in this process of 
virialization as we describe below. 

First of all, the total energy has to be conserved always, 
hence the total energy at the time of turnaround should be 
equal to that at the time of virialization. At the turnaround 
the potential energy only contributes to the total energy 
of the system whereas the point of virialization is defined 
as where the virial theorem holds i.e the kinetic energy is 
related to the potential energy T^tr = |(r^)„ir, where T 
and U are the kinetic and potential energy respectively. In 
our subsequent calculations we investigate two specific cases 
mentioned below: 

• In the first case, we assume the dark energy is inho- 
mogeneous and it virializes together with the matter. This 
assumption affects the way one calculates the turnaround 
point, as mentioned in equation (28) ( the case "Clustered 
D.E"). On the other hand, the assumption that the dark 
energy takes part in the virialization process demands that 
the virial theorem should hold for the total kinetic and po- 
tential energy of the system (i.e matter -I- dark energy). 
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Model 


Potential 


Case 


a 


b 


n 




</> 


inliom 


0.88 ±0.02 


0.34 ±0.03 


1.34 ±0.07 






horn 


0.88 ±0.02 


0.42 ±0.02 


1.25 ±0.05 






inhom 


0.855 ± 0.004 


0.383 ± 0.006 


1.10 ±0.02 


Tachyon 




horn 


0.853 ± 0.004 


0.432 ±0.005 


1.06 ±0.01 






inhom 


0.8472 ± 0.0003 


0.3867 ±0.0005 


1.001 ±0.001 






horn 


0.8462 ± 0.0001 


0.4254 ± 0.0002 


0.9824 ± 0.0005 




e'l' 


inhom 


0.850 ± 0.001 


0.387 ±0.002 


1.032 ±0.005 






hom 


0.848 ± 0.001 


0.429 ± 0.001 


1.007 ±0.003 






inhom 


0.862 ± 0.008 


0.37 ±0.01 


1.16 ±0.03 






hom 


0.859 ± 0.007 


0.43 ±0.01 


1.11 ±0.02 






inhom 


0.852 ± 0.003 


0.333 ± 0.004 


1.08 ±0.01 






hom 


0.850 ± 0.002 


0.431 ± 0.003 


1.024 ±0.005 


Scalar 




inhom 


0.8460 ± 0.0002 


0.3865 ± 0.0003 


0.9867 ±0.0007 






hom 


0.8453 ± 0.0006 


0.4241 ± 0.0008 


0.972 ± 0.002 




e'l' 


inhom 


0.8480 ± 0.0007 


0.387 ±0.001 


1.011 ±0.002 






hom 


0.8469 ± 0.0002 


0.4275 ± 0.0003 


0.9905 ± 0.0006 




PNGB 


inhom 


0.858 ± 0.006 


0.381 ±0.008 


1.12 ±0.02 






hom 


0.855 ± 0.005 


0.434 ± 0.007 


1.08 ±0.01 



Table 1. Values of the fitting parameters a, b and c by fitting equation (44) with the exact results given by equation (40) for scalar and 
tachyon models. We have also shown the errors for different parameters at 95% confidence level 



Thus, applying the energy conservation together with this, 
we get([Basilakos & Vo glisi.2007. ) 

1, 



-Ucc + Utfic 



Tc,f + UgcJ + U^aJ ~ UGc,t + U^c,i 



(34) 
(35) 



where Tc is the kinetic energy and Ugc = -3GM'^/5R 
is the potential energy for matter, M being the matter 
mass contained in the spherical over density. U,j,c = — (1 + 
3w)p0^ ^2^M jj2 ^j^g potential energy associated with the 
dark energy inside the spherical over density. Here the sub- 
scripts "f" and "t" indicate the virialization and turn around 
time respectively and "c" for inside the cluster. Using the 
above formulation, one can obtain a cubic equation which 
relates the ratio between the final (virial) Rf and the turn- 



around radius Rt, defined as the collapse factor A — 
2niA^ - (2 + n2)A + 1 = 



where 



ni = — (3it)(a/) + 1) 



i^0o/("/) 



and 



n2 



-{3w{at) 



AO -3 



(36) 



(37) 



(38) 



with pmc.t being the matter density inside the sphere at 
the turn around time while pm.t is the background matter 
density at the same epoch. 

• In the second case, we assume that the dark energy is 
homogeneous and it does not virialize inside the cluster. In 
this case, to get the turnaround point, we use eqn(28) (with 
the choice "Homogeneous DE"). On the other hand, as the 
dark energy does not virialze, its only effect is to contribute 
to the potential ener gy of the system. So ap plying the energy 
get (|Maor fc Lahavll2005l ) 



conservation, we 



t/^c + Ugc + 



RdUcc 



2 dR 



= [Ugc + U^c]t 



(39) 



(40) 



Assuming that at the collapse point, the system has viri- 
alised fully, the density contrast at virialization as a function 
of Zcou and Qmo is given by 

^ ^ /Omc,f ^ C_ (!^Y 
Pm,f A^ Vat/ 

where the relation between at and at can be estimated from 
the equation, 

^ = (41) 
da H{a)a ' ^ ^ 

together with the condition that the time needed to collapse 
is twice the turn-around time, t{ = 2tt, i.e. 



-da = 2 



-da. 



(42) 



H(a)a Jy H{a)a 

In case of a A cosmology we get an analytical solution: 

sinh-'(a/^/^^) = 2sinh-'(at^/'0^), (43) 

where vo = (1 — Q.mo) /^mo- The ratio between the scale 
factors converges to the Einstein de Sitter value (^) ~ 2^'''^ 
at high redshifts. 

In figure 5 and 6, we show the behavior of matter density 
contrast at virialization A„ir for different thawing models as 
a function of collapsed redshift Zcoii- This parameter plays 
a crucial role as an observational tool which can be directly 
applied to Press-Schechter theory for comoving number den- 
sity of virialized objects. The behavior is consistent with the 
fact that with ZcoU, ^vir decreases and settles to the value 
IStt^ for Einstein de-Sitter universe. 

Once virialized, the density contrast inside the cluster is 
slightly higher in the homogeneous case compared to the in- 
homogeneous case. This is consistent with the result at the 
turnaround where also similar things happen. This means 
that in the homogeneous case, where only matter virializes, 
the density inside the clustered objects are hi gher. This is 
in agreement with the result obtained by Maor (|Maoij|2006l ) 
for dark energy with constant equation of state. There are 
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Figure 7. Plot of the fitting function for the density contrast at virialization Ayir vs the collpase redshift (^con) for different models. 
Prom top to bottom represents tachyon with linear potential, scalar field with linear potential and scalar field with PNGB potential. For 
each case, the smooth line is for the fitting equation given by eqn (44) with the parameters mentioned in table 1, whereas the dots are 
generated by the exact numerical result for A„ir calculated for each model. 



also few interesting results. The difference in A^ir for a par- 
ticular thawing model with the corresponding ACDM model 
is larger for scalar and tachyon field with linear potential as 
well as for scalar field with PNGB potential. Also tachyon 
field with linear potential has the largest deviation from 
ACDM. This predicts that these models are easier to distin- 
guish from ACDM model by observing abundances of bound 
objects. For scalar field model, the diference between homo- 
geneous and inhomogeneous model is highest for the poten- 
tial V{(l)) = It is interesting to note that this potential 
also can act as a tracker or freezing model. For tachyon 
model, these two cases are hardly distinguishable for all of 
the potentials. 

Given the fact that the mass of the dark energy has to 
be extremely small, it is most likely that the field responsible 
for dark energy has to be homogeneous and should not clus- 
ter inside the virialized objects. If this is the scenario which 
is physically relevant, our results show that scalar as well as 
tachyon field dark energy models can be distinguished from 
the ACDM model by observing the abundances of bound ob- 
jects since the virialized density contrast for homogeneous 
cases deviate sufficiently from the ACDM model. 

We should mention, that the final results for A^ir, de- 
pend on both the background evolution as well as the fact 
that dark energy may cluster. The fact that the A^ir is 
slightly smaller for inhomogeneous case than the homoge- 
neous one for all potentials, shows that the inhomogeneous 
dark energy acts against the matter clustering. But as the 
model is quite nonlinear and involves rigorous numerical 
computations, it is very difficult to predict the percentage of 
effects coming from the background evolution as well as from 
dark energy clustering which may be an useful information 
for actual N-body simulations. 



Next, we give a fitting formula for A^ir as a function of 
collapsed redshift Zcou which is given by 

A^iriz^M) = IStt^ {a + b* eiz^oiiD (44) 

where O(zcoh) = o — r r — 1. In Table 1, we show 

the values of the fitting parameters a, b and n for different 
models assuming QmO = 0.25. We also quote the the corre- 
sponding errors for these parameters at the 95% confidence 
level. 

In figure 7, we show accuracy of our fitting function 
for thawing models with linear as well as PNGB potential 
potential. For other potentials it also works equally well. 

Knowing that models with linear potential (for both 
tachyon and ordinary scalar fields) as well as models with 
PNGB type potential deviate more from ACDM model in 
comparison to other potentials, we study the size of the col- 
lapsed objects for these potentials for homongeneous as well 
as for inhomogeneous cases. In figure 8, we show the evolu- 
tion of the perturbation collapsing radii at the present time 
{zcoii = 0) for these potentials. The radius has been normal- 
ized to the turnaround radius Rt- It has been plotted with 
respect to the scale factor normalized to the turnaround 
scale factor for ACDM model. The collapse to singularity 
is avoided assuming that after the virialization, the size be- 
comes constant. From the figure 8, for homogeneous case, 
virialization occurs earlier than ACDM model whereas for 
inhomogeneous case, it happens later. This means, homoge- 
neous models are more suitable for producing older bound 
objcets than ACDM. This gives a slight edge to homoge- 
neous dark energy models over the ACDM model. 

Although, in figure 8, the virialized radii for different 
models are practically indistinguishable but still there are 
slight differences. First of all, the virialized radii for differ- 
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Figure 8. Variation of the radius R normalized to the turnaround radius Rt w.r.t. the scale factor normalized to the turnaround scale 
factor for the ACDM model. Solid line is for ACDM whereas the dashed, dash-dotted and dotted are for scalar field woth linear potential, 
tachyon with linear potential and scalar field with PNGB potential. 



ent scalar and tachyon models are smaller than the ACDM 
model, which in turn confirms that objects are more denser 
in these models than ACDM model. This is consistent with 
the results shown in figure 5 and 6. Secondly between in- 
dividual models, virialised radius for tachyon model with 
linear potential is smaller than the other two models, e.g, 
scalar field with linear potential and scalar field with PNGB 
potential. This is also consistent with behaviour of virial- 
ized density contrast shown in figure 7. Also the virialized 
radii in inhomogeous models are slightly larger than the ho- 
mogneous counterpart, making the homogeneous one denser 
than inhomogeneous one. This is also consistent with what 
we show in figure 5 and 6. 



5 CONCLUSION 

In this work, we study the evolution of spherical overdensi- 
ties for a general class of thawing dark energy models. We 
consider both ordinary scalar fields as well as tachyon fields 
having noncanonical kinetic energy. We consider variety of 
potentials for the fields. In studying the evolution of spheri- 
cal over densities, we consider the dark energy to be homoge- 
neous as well as inhomogeneous. For the homogeneous case, 
we assume only matter virializes inside the cluster whereas 
for the inhomogeneous case, we assume both matter and 
dark energy virializes inside the cluster. 

Our main motivation is to see whether one can dis- 
tinguish thawing dark energy models from ACDM model 
studying the evolution of spherical over density. Our results 
show that although all the models become indistinguishable 
for objects collapsing earlier, for objects collapsing around 
present time, some of the thawing models deviate signifi- 
cantly from the ACDM model. To be specific, models with 



linear potential both for standard scalar fields and tachyon, 
together with scalar field with PNGB type potential have 
significant deviation from ACDM model. We have shown 
that the size of the virialised object, the density contrast at 
turnaround as well as at virialization differ significantly for 
these models from ACDM mo del. This i s consistent with the 
results obtained by Sen et al. (|Senll2O10l ) for thawing models 
considering the background evolution. Also the deviations 
are enhanced in the homogeneous dark energy case where 
matter only virializes inside the cluster, thereby making it 
more probable to be distinguished from ACDM. Given the 
fact that the mass of the dark energy field has to be ex- 
tremely small, it is safe to assume that the dark energy re- 
mains homogeneous for most of the relevant scales. Hence 
the evolution of spherical overdensities in large scale struc- 
ture formation can be a useful tool to distinguigh model 
from ACDM. 

We have also derived the fitting formula for Avir as 
function of collapsed redshift for different models and tabu- 
lated the values of the fitting parameters for different mod- 
els assuming the dark energy to be homogeneous as well 
as inhomogeneous. Although we have not shown explicitly, 
if one varies fimo, deviations are higher for values smaller 
than fimo ~ 0.25 and vice versa. This is because smaller 
the value of QmO, higher the contribution from dark energy 
at late times and hence one needs more matter inside the 
cluster for objects collapsing at late times. 

An important conclusion one can draw from this study 
is that thawing models with linear potential can have signif- 
icant deviation from ACDM model. Now, scalar fields with 
linear potential ha ve received little at tention as a viable dark 
energy candidate l|Kratochvil 1 12004| ). But it certainly de- 
serves more attention, especially for the homogeneous case. 

FinalUy, we want to stress that we study a toy model for 
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nonlinear evolution of over density where we use simplified 
assumption of spherical over-densities. In order to compare 
the theoretical results with observational data, one has to 
give away this simplified assumption and has to perform full 
N-body simulations to calculate accurately the dark energy 
signatures in nonlinear structure formation. But even with 
this simplified assumption, we find some interesting features 
in the nonlinear evolution of matter over-densities in the 
thawing class of dark energy models. This may be helpful 
for more detailed N-body simulations. 
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